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Oscillations from high energy photons into light pseudoscalar particles in an external magnetic 
field are expected to occur in some extensions of the standard model. It is usually assumed that 
those axionlike particles (ALPs) could produce a drop in the energy spectra of gamma ray sources 
and possibly decrease the opacity of the Universe for TeV gamma rays. We show here that these 
assumptions are in fact based on an average behavior that cannot happen in real observations of 
single sources. We propose a new method to search for photon- ALP oscillations, taking advantage of 
the fact that a single observation would deviate from the average expectation. Our method is based 
on the search for irregularities in the energy spectra of gamma ray sources. We predict features that 
are unlikely to be produced by known astrophysical processes and a new signature of ALPs that is 
easily falsifiable. 



PACS numbers: 14.80.Va, 98.70.Vc 



I. INTRODUCTION 

Light pseudoscalar particles appear in many exten- 
sions of the standard model. The most typical example is 
the axion, which was introduced as a consequence of the 
Peccei-Quinn mechanism to solve the puzzle of the ab- 
sence of CP violation in quantum chromo dynamics [ll, [2l ■ 
The axion is a hypothetical light particle that has a two- 
photon vertex described by the interaction term 
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where g is the axion-photon coupling constant, F is the 
electromagnetic tensor, F its dual, E the electric field, 
B the magnetic field and a the axion field. This term 
implies the possibility of photon-axion oscillations in an 
external magnetic field 0, 3 ■ This coupling is used expe- 
rimentally to search for axions that would be thermally 
produced in the Sun or axion dark matter @. In the 
case of the Peccei-Quinn axion, the photon-axion cou- 
pling is predicted to scale with the axion mass ; however, 
other models predict light pseudoscalar particles with the 
same coupling to the electromagnetic field but a priori 
unrelated to their mass Those are called axionlike 
particles (ALPs), the phenomenology of which is simi- 
lar to standard axions. Astrophysical environments can 
offer ideal conditions for photon-ALP oscillations, with 
the possibility of long baseline experiments involving ma- 
gnetic fields Q. Progress over the last decade in 7-ray 
astronomy allowed one to consider searching for the im- 
prints of 7- ALP oscillations in the energy spectra of high 
energy 7-ray sources @ . The effect of 7- ALP oscillation 
is usually assumed to be twofold : it is expected to induce 
a dimming of the fiuxcs above a given threshold [l^, [HI , 
and possibly decrease the gamma-ray pair production re- 
lated opacity at high energy. The opacity can be that of 
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the intergalactic medium [l2|, |l3| or within the sources 
themselves [3] . A crucial point is the turbulent nature of 
the magnetic fields the photon beam travels through. It 
implies a consequential randomness in the prediction of 
the observable effects. This has been pointed out in [l5| 
in the case of the change of opacity due to 7- ALP oscil- 
lations. The authors of [l^, [la| showed that because of 
the random nature of the intergalactic magnetic fields, 
the effect of 7- ALP mixing should be very different from 
one source to another. Such an observable is then useless 
to perform ALP searches through the observation of a 
single source, leaving only the possibility of a population 
study in order to average the effect over many sources. 
This type of study has been conducted in e.g. sho- 
wing a hint for an anomaly in the transparency of the 
Universe. Though rapidly increasing with the advent of 
the last generation Cherenkov telescope arrays such as 
HESS, MAGIC and VERITAS, there arc only a handful 
of high energy sources that are effectively concerned with 
extragalactic absorption. It is thus interesting to point 
out some effect of the 7- ALP mixing that does not rely 
either on stacking or averaging, in order to exploit obser- 
vations of single sources. Here for the first time an effect 
is pointed out that potentially applies to single observa- 
tions. This article is organized as follows. First, we briefiy 
recall the formalism of 7- ALP mixing and apply the re- 
sults to a single coherent magnetic domain. As a second 
step we show the results of a simulation of photons trave- 
ling through a set of magnetic domains. In particular we 
show that contrary to what is stated in the literature, a 
sharp drop in the energy spectrum of high-energy 7-ray 
sources is not a robust observable and is not what should 
be searched for. Actually the 7- ALP mixing would pro- 
duce an anomalous dispersion of the spectra, which would 
no longer appear as smooth in a limited energy range. 
We then give an explicit example of how the effect could 
appear in the data, in the case of a specific situation, na- 
mely an extragalactic TeV emitter whose photons travel 
through the intergalactic magnetic field (IGMF), and we 
discuss the robustness of the method. 
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II. THE PHOTON/ AXION SYSTEM IN A 
MAGNETIC FIELD 



The 7- ALP system is described following the approach 
of Q ■ A three-state wave function is used with two states 
of polarization for the photon and one for the ALP. Let 9 
be the angle between the magnetic field direction and the 
photon momentum. Since only the B component trans- 
verse to the propagation couples photons and ALPs, the 
strength of the magnetic field involved in the coupling is 
B sinO. Moreover, for parity issues, only one polarization 
state parallel to the field is involved in the interaction. 
This is accounted for by introducing the angle (j) between 
the transverse component of the field and the direction of 
the polarization sate Ai . The 7- ALP system is then pro- 
pagated using the following linearized equations of mo- 
tion assuming relativistic axions : 
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with the mixing matrix 

Aii-iAabs A12 AbC0S( 
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where Ab = gB sin 6/2 is the coupling term, and Aa = 
—rril/2E is the ALP mass term. Here we neglect the Fa- 
raday effect and the vacuum Cotton-Mouton term, as 
the low magnetic field strength considered in the follo- 
wing makes the corresponding contribution irrelevant for 
this study. This implies A12 = A21 ~ and that the 
other diagonal terms are An = A22 = — ajpj/2_E, ujp\ 
being the plasma frequency accounting for the effective 
photon mass. As in |18| . absorption of photons on their 
way is introduced with the Aabs = t/2s term where r is 
the optical depth assuming a propagation over a domain 
of size s within which the opacity is homogeneous. Be- 
cause of that term, the matrix is no longer Hermitian and 
unitarity is lost. In the following, this term will be used 
to model the absorption of photons on the extragalactic 
background light (EBL) while propagating in IGMFs. Af- 
ter diagonalization of the mixing matrix, the equations 
of motion can be analytically solved and the transfer ma- 
trix of the system is obtained. The probability of 7 — a 
conversion after crossing one coherent magnetic field do- 
main of size s in the simplest case, without absorption 
and neglecting the plasma term, yields 
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with Aqsc = \/ Aa + 4Ag. The energy dependence of the 
mass terms in Aqsc implies an energy threshold above 
which the conversion becomes efficient. 
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Figure 1: Survival probability of an unpolarized photon as a 
function of the energy for three values of S. 



nies being the effective ALP mass in the presence of 
charges (e.g. in a plasma). For E <C £'thr, Aqsc ^ Ab and 
then no conversion occurs. For E ^ Ethr spectral oscilla- 
tions happen due to the energy dependent sin^ Aoscs/2 
term. For E ^ Ethr , Aqsc ~ Ab and the conversion pro- 
bability is no longer energy dependent. The conversion 
probability of Eq. 2] can be parameterized in terms of 
Ethr and 



6 ^ gBs sine /2 
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instead of Bs'mO and s. sm^S/2 is then the conversion 
probability at very high energy (VHE, E ^ -Ethr)- The 
condition required for a significant conversion to occur, 
(5 > 1, is similar to the Hillas criterion for the acceleration 
of ultra high energy cosmic rays, as pointed out in [loj . 
Figure [1] shows the evolution of the photon survival pro- 
bability as function of the energy for three different va- 
lues of S. For allowed large IGMF values of order 1 nG, an 
ALP mass of 2 neV and a coupling g — 8 x 10~^^ GeV~^ 
at the limit of current experimental constraints [H , i^thr 
lies at about 1 TeV. The asymptotical value of 1 — P^^a 
gives the level of dimming of the photon fiux (indepen- 
dently of eventual additional EBL absorption). One can 
see that this attenuation is hardly predictable given the 
uncertainties on the environmental parameters, B and s, 
as it depends sinusoidally on the value of S. 

In astrophysical environments, magnetic fields are 
usually not coherent. In the case of a propagation through 
a turbulent magnetic field, the beam path can be divided 
into coherent domains of size of the coherence length of 
the field (the validity of this simple model is discussed in 
Sec. lIV[) . For each domain, a transfer matrix is generated 
with a random orientation of the magnetic field yielding 
a specific value of 5. The total transfer matrix associated 
with this realization of the turbulent magnetic field is the 
product of all individual transfer matrices. The spectral 
shape of the global conversion probability for one single 
realization is the result of the interference of all oscillation 
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patterns such as those displayed in Fig.[TJ As the pseudo- 
period is different in each domain, the photon survival 
probability has a very complex energy dependence. As 
an illustration, the survival probability of a photon from 
a source at redshift z = 0.1 traveling through a single 
realization of a 1 nG IGMF with coherent domains of 
size So = 1 Mpc is displayed in Fig. [21 A plasma density 
of He = 10~^cm~'^ typical of the intergalactic medium 
is assumed. In this condition and for ALP masses of or- 
der neV, TOoff = nia.- For illustration, the upper panel 
shows the survival probability without absorption on the 
EBL, whereas the lower panel results include this effect. 
Conservatively, the EBL density model used here is the 
lower limit model from [l^. To account for redshifting, 
a flat ACDM Universe with {Ha, f^m) = (0.73, 0.27) and 
Hq = 71 km/s/Mpc is assumed. Here the dashed red line 
is the prediction without ALPs, so that the dimming is 
only due to EBL. From Fig. [2] one can see that the pre- 
diction of the model including ALPs is the presence of 
a significant level of noise in the energy spectrum over 
one decade or so around E'thr- Because of the unknown 
nature of the orientation of the magnetic field within the 
domains, the exact shape of the spectrum in this region is 
unpredictable. However, as we shall see in the following, 
the noise level is a prediction of the model. This predic- 
tion significantly differs from what usually appears in the 
literature, namely a smooth transition between no dim- 
ming below Ethi and a fixed level of attenuation above it. 
It has been shown in [2^ that the averaging over a large 
number of realizations of N domains in each of which the 
conversion probability is Pq yields an overall conversion 
probability 

^7- = J(l-e"^^^°) . (7) 

This means that the effect as it has been studied so far 
is valid for an average over a collection of sources. In the 
case of the observation of one source only, if N is very 
large and the energy spectrum is binned, then the smooth 
behavior can be retrieved in principle. In practice N is 
not large enough, as we shall see in the following. The 
results presented in Fig. [2] are obtained with a single rea- 
lization. By averaging the results of Fig. [2] over a large 
number of realizations, the value given by Eq. [7]is retrie- 
ved. From one realization to another, only the orienta- 
tions of the magnetic fields vary ; the number of domain 
and their sizes are kept fixed. 

Note that above 5 TeV the survival probability for this 
specific realization is higher with ALPs than with EBL 
only. This is the so-called opacity effect, because pho- 
tons are untouched by the EBL as they travel disguised 
as axions, the Universe appears to be more transparent. 
This result should be taken with care, however because, 
as shown in [31, there exist realizations of the IGMF 
where the opposite effect is obtained, basically when most 
ALPs do not convert back to photons before detection. 
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Figure 2: Photon survival probability as function of the 
energy for a realization of a source at z — 0.1 using B = 
1 nG, So = 1 Mpc, 5 = 8 X 10"" GeV"^ and nia = 2 neV 
without absorption (upper panel) and with EBL absorption 
(lower panel). 



III. OBSERVATIONAL EFFECTS 

The experimental relevance of the proposed signature 
is now studied in the particular case of a source at red- 
shift z = 0.1 for the same parameters as above. The in- 
trinsic spectrum of the source is simulated following a 
log-parabola shape with an integrated flux in the TeV 
band at the Crab level. A 50 h observation time is assu- 
med with an energy resolution of 15 % and assuming a 
constant effective area of 10^ m^, these values being ty- 
pical of current generation Cherenkov observatories. The 
intrinsic spectrum is convolved by one randomly gene- 
rated photon survival probability and eventually binned 
to obtain the spectrum that would be observed in this 
model. The result of this simulation is displayed on the 
left panel of Fig. [31 A fit of the simulated experimen- 
tal data by a log-parabola shape convolved with EBL 
absorption is also shown, as it would be performed by 
observers. In the right panel of Fig. [3|are displayed the 
residuals of that fit. It appears that in the case without 
ALPs the residuals would evenly spread around whe- 
reas these residuals would show anomalously strong and 
chaotic deviations from in the case of 7-ALP mixing. 
This is the expected signature of ALPs in the spectrum, 
induced by the noisy spectral shape of the photon survi- 
val probability. 

The approach considered here corresponds to what an 
observer would do. First one would fit a smooth shape 
(be it a log-parabola, a broken power law or an emission 
model inspired one) and then pick up the shape provi- 
ding the best and a decent residual distribution {e.g. 
centered on zero, without obvious biases, etc.). The cru- 
cial point is that the observer would fit the data with a 
smooth function. This is motivated by the fact that no 
TeV source emission model predicts spectra with local 
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Figure 3: Simulation of the observation of one 7-ray source at 2=0.1, with the effect of 7-ALP mixing (left), and distribution 
of the residuals of fits to a conventional model and a model with ALPs (right). 



extrema or a noisy shape. So this result holds for any 
smooth spectral shape provided it gives the best pos- 
sible fit. For a given observation, a given emission mo- 
del and a given fit function, the noisy energy range and 
the variance of the residuals is a prediction of the ALP 
model. This is illustrated in Tab. HI virhich displays the 
variance of the fit residuals under different hypotheses : 
no ALP and two values of g, with an ALP mass still 
yielding i^thr = 1 TcV. The exact value of the variance 
of the fit residuals depends on the analysis that would 
be performed, in particular the energy range chosen by 
the observer. The use of the variance of the fit residuals 
is only an example, as observers might choose to use a 
more sophisticated estimator of the noise in theirs data 
sets. Additionally, because of the random nature of the 
predicted effect, it is important here to verify that the 
scatter in the prediction among realizations is smaller 
than the effect itself. The predicted uncertainty on the 
variance of the fit residuals due to the random nature 
of the prediction is also shown in Tab. |T] for each consi- 
dered scenario. These results are obtained by averaging 
over 5000 realizations. It appears that for the considered 
parameters, the effect is significant, as one can see that 
the variance anomaly in the presence of ALP is predicted 
to be significantly above the conventional value. The ob- 
servation of a Crab-level source for 50 h was chosen here 
as an illustration. Actually one finds that for the same 
redshift and energy range, the effect would still be visible 
but with less significance by observing only 5 h. 

The observational signature that is discussed here oc- 
curs for energies around -Ethr given in Eq. [5] Therefore 
the range of accessible ALP parameters with this method 
depends on the value of the magnetic field and the energy 
range of the experiment. For instance considering TeV 7- 
rays and nG IGMF, the above results show that a typical 
lACT would be sensitive to ALPs with g - 10"" GcV^^ 
in a range from 0.1 neV to 10 neV. In that range of 
mass the most stringent constraint currently comes from 



Model Variance of the fit residuals 

No ALP 0.04 ±0.01 

g = 10~", m = 0.7 0.11 ± 0.04 

g = 8 X 10~", m = 2 0.20 ± 0.05 

Table I: Values of the RMS of the fit residuals to mock data 
with different assumptions for g and m (in units of GeV~^ 
and neV resp.), for constant size magnetic field domains. 



the CAST helioscope with an upper limit on g of order 
10~^° GeV~^ 0. So in principle this method should al- 
low improving current constraints in this range of mass. 
To go to larger masses, one has to consider larger magne- 
tic fields (in principle the method discussed in this article 
is valid for any 7-ray source behind a turbulent magnetic 
field) and/or higher energies, as the relevant mass for a 
given g goes as \/E x B. 

IV. DISCUSSION 

One important point is that should anomalous disper- 
sion be observed some day, one would know how to fal- 
sify the interpretation in terms of new physics. This can 
be done for instance by observing the same object with 
more exposure. If the ALP interpretation is wrong, lo- 
cal extrema would not hold and all the residual points 
would be redistributed around zero. If the interpretation 
is correct though, two effects arc predicted due to the in- 
creased statistics : i) the significance of the deviant bins 
would strengthen, and ii) irregularity would disappear 
at VHE as expected from ALP models. The first point is 
justified by the fact that a magnetic field that is coherent 
over a scale s should remain coherent over times of oder 
s/c. For scales of order 1 Mpc as relevant here, this time 
scale is of order 3 x 10^ yrs. Concerning possible effects 
that could produce similar irregular spectra, one could 



5 



Model Variance of the fit residuals 

g = IQ-", m = 0.35 0.18 ±0.05 

ff = 8 X 10~", m = 1 0.42 ±0.14 

Table II; Values of the RMS of the fit residuals to mock data 
with different assumptions for g and m (in units of GeV^^ 
and neV resp.), in the case of a Kolmogorov-like turbulent 
magnetic field. 

imagine a complex landscape of background UV-IR pho- 
tons that produces non-trivial absorption features in the 
energy spectra and mimics the effect. In the event of a 
positive detection, this would therefore require studying 
the effect over more sources and how it depends on z for 
instance. For observers interested in putting constraints 
on ALP models, though, this is not an issue since such an 
effect would add up to the irregularity of the spectrum 
and by no way it could cancel it. 

The modeling of the IGMF as it is done here with do- 
mains of same sizes is the simplest model one could think 
of. It has been used here as it is widely used in the li- 
terature. To describe more precisely the magnetic field 
turbulence, it is possible to account for the power distri- 
bution of the modes. The turbulent field can be modeled 
as a Gaussian random field with each Fourier mode pro- 
portional to some power of the wave number In the 
generic case of isotropic and homogeneous Kolmogorov- 
like turbulence, a = 5/3. As shown in 0], this leads to a 
variation of the rms intensity of the magnetic field B as 
a function of the scale s such that B oc s^/'^ . 

Before discussing the effect of such a magnetic field 
on 7-ray source spectra, let us remark that magnetic 
fields that are coherent on small scales should have ne- 
gligible effects on the spectrum in comparison with the 
larger scales. For small values of (5, Poj the conversion 
probability over a scale s, is expected to be of order 
(5^/2 ^ g^B'^s^ 1% (see Eq. [5]). In that case, the averaged 
formula of Eq. [7] reduces to P-^^a — N5'^ /2. All in all, for 
a given <?, this probability is proportional to NB'^s^ . If Ps 
is the probability of photon conversion for modes of size 
s, given the above mentioned law for the magnetic field 
strength, one gets Ps/io ~ 2.5% x P^ for the conversion 
probability in a magnetic field mode corresponding to a 
scale s/10. This means that the small scales rapidly be- 
come irrelevant for this study and one can safely consider 



that the largest scales contribute the most in the power 
distribution of modes. Concerning larger scales, the ef- 
fect on the noise level is limited by the ratio between the 
considered scale and the distance to the source. Speaking 
in terms of domains, if there are only a few equivalent do- 
mains, little interference will happen and then the noise 
in the energy spectra will have wider fluctuations. 

To be more quantitative, the study of Sec. IIIII has 
been repeated using a Kolmogorov-like turbulent magne- 
tic field inspired by the modeling used in [2l| . As for the 
previous study, 5000 realizations of turbulent magnetic 
field are performed, with wave numbers ranging from 0.1 
Mpc to 100 Mpc, and a rms intensity of i? of 1 nG at 100 
Mpc. The exact same kind of noise in the 7-ray spectra 
is obtained. To illustrate this, the results of these simula- 
tions are shown in Table |lT]; in particular, the variance of 
the fit residuals is still larger than in the no- ALP situa- 
tion, in a statistically significant way. It has been checked 
that this results is stable when larger scales are used for 
the lowest wave number. 



V. CONCLUSION 

In this study we showed a new possible signature of 
7- ALP mixing in the form of an anomalous dispersion in 
the energy spectra of 7-ray sources. The smooth-noisy- 
smooth alternation behavior in the energy spectrum is 
a peculiar prediction of ALP models that could hardly 
be mimicked by known astrophysical processes. It has 
been shown that this effect can be used to constrain ALP 
models from the observation of single sources. An expli- 
cit example has been given in the case of oscillations in 
IGMF ; however, such a signature can be searched in any 
source for which a turbulent magnetic field is present 
along the line of sight. 
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